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OBJECTIVES 


• Develop an efficient and robust algorithm for multi-phase 
chemically reacting flows at all speeds, with emphasis on 
low Mach number flows. 

• Calculate turbulent spray combustion flow in a gas tur- 
bine combustor. 


MOTIVATION 


Many reacting flows in propulsion devices cannot be efficiently calculated by 
modern compressible flow CFD algorithms, e.g., 

- rocket motor — wide range of Mach numbers, from near zero 
velocity at closed end to supersonic at nozzle exit. 

- gas turbine combustor — low subsonic velocity, but large density 
variation precludes incompressible approach. 

Most low-speed reacting flow codes based on TEACH-type technologies - 
inefficient and lack of robustness for complex flows. 

Tremendous progress made in high-speed compressible flow CFD in past 
two decades. Extending application range to low-speed flow regime highly 
desirable. 


OUTLINES 


• GOVERNING EQUATIONS 

- Gas-Phase Equations 
— Liquid-Phase Equations 
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• FUTURE PLAN FOR ALLSPD CODE 



GOVERNING EQUATIONS 


• Gas-Phase Equations 
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ALLSPD MATN FEATURES 


Present Capabilities 

• 2-D and Axisymmetric geometries 

• Second-order central difference for both inviscid and viscous terms 

• Fully coupled, fully implicit algorithm 

• Efficient convergence for wide range of Mach numbers (from M < 10 to 
supersonic) 

• Finite-rate chemistry, realistic thermophysical properties 

• Multi-block and body-fitted curvilinear coordinates for complex geometries 

• k — e turbulence model . 

• Stochastic liquid spray model (dilute spray), vortex model for droplet internal 

circulation and diffusion 

• Still a research code, require experience and knowledge of CFD and flow 
physics to use. 

Future Plans 

• Development of a more efficient solver 

• Extension to 3-D 1 

• PDF model for turbulence/chemical reaction closure 

• Thermal radiation model 

• Detailed soot and NOx kinetic models 

• Multi-grid and unstructured grid capabilities 

• Dense spray and high pressure (near-/super-critical) spray models 


Difficulties with Compressible Flow Algorithms at Low Mach 
Numbers 


• Disparities among system’s eigenvalues (stiffness), u, u + c, u - c, re- 
sulting in significant slowdown in convergence rate. 

• Singular behavior of pressure gradient term in momentum equations 
as Mach number approaches zero, 
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As Mach number is decreased, pressure variation (Ap* oc M 2 ) becomes 
of similar magnitude as roundoff error of the large pressure gradient 
term (p^-yM 2 ). 


METHOD OF APPROACH 


Pressure Singularity Problem 

• Pressure decomposed into two parts: 


p= P 0 + Pg 


Pg replaces p in momentum equations and retains p g as one of the 
unknowns. 

• Employs conservative form of governing equations, but uses primitive 
variables 

(p. Jt u,v,h,Yi) 


as unknowns. Conservation property preserved and pressure field ac 
curately resolved for all Mach numbers. 


Eigenvalue Stiffness Problem 

• Pressure rescaled so that all eigenvalues have the same order of mag- 
nitude. Physical acoustic waves removed and replaced with pseudo- 
acoustic waves which travel at speed comparable to fluid convective 

velocity. 


EIGENVALUES RESCALING 


A = {/, U, l [t/(l + ± - ^) 2 + 40(«i + “!)]. U,U--- 


<*! = {„ « 2 = £y. l/ = ai u + a 2 v. 

For well-conditioned eigenvalues, scaling factor /? taken to be 


P = u 2 + w 2 . 


LIQUID-PHASE EQUATIONS 


Droplet Motion Equations 
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• Droplet Heat and Mass Transfer Equations 
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Droplet Internal Temperature Equations 
(Vortex Model) 
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NUMERICAL ALGORITHM 


Gas-Phase - ALLSPD code 
Liquid-Phase 

— Droplet motion equations (ODE) - Runge-Kutta method. 

- Droplet internal equations (PDE) - implicit method 
(Thomas algorithm). 

— Determination of spray time step for integration. 

- Stochastic separate flow model. 

Interaction Between Two Phases 


SPRAY TIME STEP 


Droplet Velocity Relaxation time ( t r ) 
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Droplet Life Time ( t\ ) 
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Droplet Surface Temperature Constraint Time ( t s ) 
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• Local Grid Time Scale ( t g ) 

• Turbulent Eddy-Droplet Interaction Time (fj) 
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INTERACTION BETWEEN TWO PHASES 


1. Initialize gas and liquid phase variables. 

2. Solve liquid-phase equations. 

3. Evaluate spray source term, Hj. 

4. Solve gas-phase equations and update gas-phase 
variables. 

5. Update spray source term, H/ ? 

No, go to step 4. 

Yes, go to step 2. 


NUMERICAL TEST RESULTS 


Turbulent Backward-Facing Step Flow - non-reacting. 
Evaporating Turbulent Spray Flow. 

Gas Turbine Spray Combustion Flow. 
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Convergence history tor both non-combustion and spray combustion cases for the gas 
turbine spray combustion flow. 
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CONCLUSION 


• ALLSPD code efficiently coupled with the SSF spray 
model. 

• Satisfactory convergence property for flows with and 
without spray. 
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